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Increasing demands for electric field-tunable electric, 
magnetic, and orbital (EMO) materials has renewed inter- 
ests in ferroelectricity and its coupling to EMO properties. 
Materials where a spontaneous electric polarization P, 
which naturally couples to and is switched by an electric 
field, strongly interacts with EMO degrees of freedom pro- 
vide a platform to realize novel electric- field-controllable, 
low-power multifunctional devices such as ultra-fast Mott- 
based devices.G^ In the versatile class of complex ABO^ 
perovskite oxides (Fig. [T]) , ferroelectric polarizations are 
usually induced by polar displacements of second-order 
Jahn-TeUer (SOJT) active cationsPtll] on the B-site (com- 
monly dp transition metals ions) and/or on the A-site 
(lone-pair active cations). The cation displacements that 
induce ferroelectricity, ho wever , are often incompatible 
with and/or weakly couplecP^I^ to EMO-d erived material 
properties such as electronic bandwidthSjISEH magnetic 
interaction JiSEII and critical transition temperatures .liSUSl 
Rotations of corner-sharing -BOg octahedra, however, di- 
rectly alter these macroscopic properties,!^ because they 
buckle the inter-octahedral B-O-B bond angles which 
mediate the interplay of EMO degrees of freedom. Thus 
in perovskites, although the vast number of chemical 
compositions available (Fig. [l^) facilitates nearly every 
conceivable material property, the local polar cation dis- 
placements necessary for ferroelectricity often occur inde- 
pendently of the property-controlling and pervasiv^^HM! 
BOq rotations in the extended framework — rotations ac- 
count for nearly 75% of all rhombohedral. Fig. [TJd, and 
orthorhombic. Fig. [l]:, structures. 

The strong coupling of oxygen octahedral rotations to 
the EMO properties in the perovskite structure presents 
an opportunity to create novel multifunctional materi- 
als that respond to external electric fields. The steric 
linkage of the -BOg framework, however, constrains ro- 
tations of octahedra in the same plane to have equal 
magnitudes but opposite rotational "sense," thereby pre- 
serving inversion symmetry.'^ Thus collective rotation 
patterns do not induce an electric polarization nor di- 
rectly couple to an electric field. To circumvent this 
reality, synthetic chemical routes have focused on turning 
the centric BOq octahedra into acentric structural units 
by partial halide substitution.'^SHll! Progress in cis- and 
frons-ordering of the halide anions throughout all acentric 
units, however, remains slow and ferroelectric switching 
even more challenging. ^SiHI j^gogntly, a new strategy in- 
volving the layering of perovskites blockj23ilIl has lead 
to the realization of ferroelectric systems whose electric 
polarizations can be completely accounted for by centric 
rotations. The mechanism has been referred to as hybrid 
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FIG. 1. The ubiquitous octahedral rotations available to 
perovskites regardless of chemistry preserve inversion symme- 
try, (a) The adaptable ABX^ perovskite structure permits a 
vast number of chemical compositions, indicated by the color- 
coded periodic table, and therefore also exhibits nearly every 
conceivable materials property. The oxide material class, es- 
pecially the compositions studied here (listed below), exhibits 
a network of corner-connected B0(, octahedra which are often 
distorted from the cubic structure shown in (a). The perva- 
sive low symmetry (centrosymmetric) rhombohedral (b) and 
orthorhombic (c) structures that result are due to rotations 
of the BO% building blocks. Following Glazer notation,'^ the 
a~a~a~ octahedral tilt pattern in (b) consists of rotations 
of adjacent BOe octahedra that are out-of-phase (-) and of 
equal magnitude along all Cartersian directions. The a~a~c'^ 
tilt pattern (c) is similar as it possesses out-of-phase rotations 
about two directions with equal magnitude (perpendicular 
to the projection illustrated), but exhibits in- phase rotations 
(4-) about the remaining Cartesian direction with different 
amplitude (page normal). Periodic table adapted after Ref. I24l 



improper ferroelectricity, however, no comprehensive de- 
sign prescription applicable to multiple chemistries has 
been proposed. 

In this work, we present a materials design strategy to 
realize ferroelectricity in layered perovskite superlattices 
constructed from two different ABO3, materials, neither 
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FIG. 2. Illustration of the Chemical Criterion for rotation-induced ferroelectricity in layered perovskite superlattices 
constructured from two different ABO3 perovskite materials. In (a) bulk ABOs perovskites inversion (I) centers are found on 
both the A- and B-sites; the highest site-symmetry operator being a three- fold rotoinversion (3). Cation ordering in layered 
perovskites, however, lifts the inversion centers on the B-site (leading to a 4- fold rotation) in the A/ A' layered perovskites (b) 
and on the yl-site in the B/B' (c) structures. Inversion remains through the ^ operation found on the remaining A-site and 
B-site, respectively. Since rotations of octahedra preserve the inversion on B-sites yet can remove it on the A-sites, only A/A' 
support this form of hybrid improper ferroelectricity. 



of which are ferroelectric in bulk. The ferroelectricity is 
of a hybrid improper fornij'SSHM! where the polarization 
P can be induced by the coexistence of two octahedral 
rotation patterns of different symmetry. The design 
rules are constructed from group-theoretical methods 
combined with principles of perovskite crystal chemistry 
and are therefore completely general. In particular for 
{AB03)i/{A' B'03)i superlattices, a layered A/A' cation 
ordering, designated as a chemical criterion, is required 
to realize this new form of ferroelectricity, whereas 
optimal materials selection, corresponding to an energetic 
criterion, is based solely on the rotation patterns present 
in the parent single phase ABO3 and A'BO^ bulk 
perovskites. We find that the parent materials do not 
need to exhibit ferroelectricity, but rather that each 
constituent should have a tendency to distort in the 
orthorhombic Pnma (a~a~c'^ in Glazer notation, see 
Fig. [ij;) octahedral tilt pattern in bulk. Using ab initio 
first-principle techniques we test the design rules against 
artificially layered gallate and aluminate perovskites. 
(Note that our focus on dP cations is to avoid the known 
challenges of first-principles theory in treating materials 
with partially filled d- and/or /-shells and is not a 
limitation of our strategy.) These materials design 
rules are capable of rapidly guiding the discovery of 
numerous unknown multifunctional materials as the 
A /A' cation layered orderin g is highly amenable to 
advanced solid-state synthesis^^l^ and layer-by-layei^ 
deposition techniques. Additionally, the targeted a~a~c+ 
octahedral tilt pattern is the most common type observed 
in perovskites. 

Structure-property guidelines 

Hybrid improper ferroelectricity (HIF) has been sho wn to 
arise from a peculiar trilinear lattice coupling ternpSlMI 
in the thermodynamic free energy, gP(QiQ2), where the 



polarization P is induced by the product of two trans- 
lation symmetry- breaking lattice modes, Qi and Q2, of 
different symmetry. The interaction strength of the cou- 
pled modes is described by the temperature-independent 
coefficient g. The principal design challenge is identifying 
the universal structural-chemical requirements within a 
class of materials that allows the non-polar lattice modes 
to induce ferroelectricity. Here, we focus on octahedral 
rotation-induced ferroelectricity in perovskite-structured 
oxides due to the abundance of Q modes describing the 
BOq connectivity. 

Our crystal-chemistry design approach is to consider 
the archetypal five-atom cubic ABO^ perovskite as a basic 
chemical unit, interleaving any two perovskites to form a 
layered bi-color superlattice structure. [It is worth remem- 
bering that neither Glazer octahedral rotation patterns 
in bulk perovskites nor cation layering of the A/ A' (or 
B/B') sites in such bi-color superlattices are capable of in- 
dividually producing noncentrosymmetric structures.] We 
consider (AB03)i/(A'B03)i and {ABO:^)^/ {AB'0^)i 
superlattices, corresponding to the bulk-like composi- 
tions {A, A')B206 and A2iB, B')Oq (Sec. SI) with layered 
cation-orderings (Fig. [2]) , as the simplest structures that 
allow us to formulate materials selection rules in terms of 
the properties of the two hulk constituents; The chemi- 
cal composition of the basic perovskite units considered, 
and the energetic or lattice dynamical properties of the 
five-atom perovskite units. 

Chemical Criterion: The synthetic perovskites require an 
A/ A' layered cation ordering for the octahedral rotations 
to induce ferroelectricity. The crystallographic rational 
for this criterion is as follows: In the ABO-^ perovskite 
structure both the A-site and the i?-site positions have 
local inversion symmetry (Fig. [2^). Rotations can lift 
the inversion center on the ^-site, but always preserve 
the centricity (inversion on the i3-site) of the octahedra. 
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FIG. 3. The primary rotational modes Qi (a a c ) and Q2 
(a'a'c") targeted to combine in the A/A' superlattices to 
induce hybrid improper ferroelectricity. The combination of 
these rotations produces the common orthorhombic (a~a~c^) 
tilt pattern found in perovskites. 



A/ A' {B/B') layered cation ordering removes the inver- 
sion center on the i?-site (A-site) as shown in Fig. 
and c. Therefore the combination of A/ A' {B/B') cation 
order and rotations allows for (forbids) hybrid improper 
ferroelectricity. 

This heuristic understanding is verified by our group 
theoretic methods. We find seven unique tilts patterns 
obtained from pairs of centric BOq octahedral rotations 
patterns in the paraelectric {A, A')B20e structure. Each 
combination leads to a P{QiQ2) term in the free energy 
(Table S5) allowing for rotation-induced ferroelectricity — 
neither SOJT cations nor polar cation displacements are 
required. Note that the form of the trilinear term indicates 
that the combined rotation pattern is effectively polar; 
once the two rotations are present, there is only one 
direction for the induced polarization, i.e., turning on 
the two rotations is analogous to turning on an electric 
field. In contrast, we find no trilinear invariants in the 
paraelectric A2{B,B')0q structure and therefore hybrid 
improper ferroelectricity is forbidden. 

The chemical criterion is insufficient to guarantee that 
the octahedral rotations induce a polarization, i.e., that 
Q1Q2 P- Rather, the form of the trilinear coupling, 
P{QiQ2), compels only the mutual coexistence of such 
lattice distortions, which are commonly incompatible with 
one another, if any two are unstable.l^ 

Energetic Criterion: The ground state structure of the 
ABO3 and the A'BOs building blocks must each contain 
the corresponding Qi and Q2 rotation patterns. This re- 
quirement provides an optimal condition for hybrid im- 
proper ferroelectricity: it ensures the two octahedral ro- 
tations in the synthetic superlattice dominate the energy 
landscape over other competing instabilities, indicating 
they drive the transition to the polar structure. Said 
another way, the ground state structure of both bulk 
chemical units should possess the rotation pattern which 
results from the specific combination of octahedral rota- 
tions modes, Qi and Q2, that are desired to be present 
in the 1/1 ABOs/A'BOs superlattice. 

This criterion stems from the fact that the rotation 



patterns are fully coherent in the synthetic A/ A' 
structures. A rotational instability that appears in 
only one of the constituents, although likely capable 
of inducing the same rotation pattern in the second 
chemical unit would produce a tilt pattern with overall 
smaller octahedral rotations and modulated rotation 
angle amplitudes. Such compound tilt patterns are 
detrimental to cooperative ordering of local dipoles and 
could hinder octahedral rotation-induced ferroelectricity. 

First-principles guided materials design 

Here we apply our criteria in order to select bulk single 
phase perovskite oxides with suitable chemical compo- 
sitions and rotational patterns for integration into the 
1/1 ordered ABO3/ A' BO3 perovskites. The first design 
criterion is necessarily satisfied by restricting our investi- 
gation to this composition. Next we apply group-theoretic 
methods and consider all symmetry-adapted normal mode 
patterns (Table S5) that produce octahedral tilt patterns 
in the paraelectric {AB03)i/{A'B03)i structure. This 
analysis reveals many different possible invariants (Eq. 
SI) that permit the octahedral rotations to induce ferro- 
electric polarizations through the trilinear PQ1Q2 term 
in the free energy. To narrow the search, we focus on the 
lattice modes Qi = a"a°c+ and Q2 — a^a^c'^ (Fig. |3|. 
These rotational modes are prime distortions to target 
because they are ubiquitous in single phase perovskites; 
they produce orthorhombic Pnma perovskites with the 
a~a~c'^ octahedral tilt pattern. 

We choose to focus on five different gallate and alu- 
minate perovskites: YAIO3, LaAlOa, BiAlOa, LaGaOa, 
and YGaOs. We evaluate the lattice stability from cal- 
culations of the force constants throughout the Brillouin 
zone of these bulk compounds (Fig. S2) and perform (zero 
kelvin) full structural relaxations in various isotropy sub- 
groups corresponding to the condensation of one or more 
of the unstable modes at the high-symmetry points within 
the Brillouin zone. In Table |T] we enumerate the most 
unstable phonon modes at the high symmetry points in 
the tetragonal Brillouin zone of the non-polar (PA/mmm) 
bulk phases and identify the rotational pattern of the per- 
ovskite ground state. Note, although some of these bulk 
compounds show polar T-point instabilities,'^ none of 
them exhibit ferroelectric perovskite ground state struc- 
tures in the bulk except for i?3c BiAlOa. 

The combination of our targeted octahedral rotation 
modes, Qi — a°a"c+ (Af-point) and Q2 = a^a^c° (A- 
point), determine the ground state tilt patterns of three 
of the five compounds: LaGaOa, YGaOs and YAIO3 
all exhibit the orthorhombic a^a^c^ rotation pattern at 
room temperature (Fig. ij:). In contrast, LaAlOs and 
BiAlOa display the rhombohedral a^a^a^ tilt pattern. 

We consider four representative 1/1 superlattice test 
cases based on these bulk candidate perovskites (Table [TT|: 

(1) YAlOa/YGaOa.— Both YAIO3 and YGaOa exhibit 
strong a^a'^c^ and a~a~c^ instabilities, i.e., the 
targeted orthorhombic rotation pattern, but the 
B/B' ordering does not satisfy the chemical crite- 
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TABLE I. Phonon modes of the single phase paraelectric alu- 
minates and gallates used in the selection of bulk perovskites 
to combine and realize rotation-induced ferroelectricity. Cal- 
culated frequency and distortion-types of the most unstable 
phonon modes for the reference PA/mmm bulk phases of the 
constituent superlattice materials. Imaginary frequencies indi- 



cate energy loweriuj 


J instabilities. 


The high- 


-symmetry fc-points 


have the following 


wave vectors: 


r = (0,0,0), Af=(iiO), 


i?=(0,|,|), andA=(i,|,|). 






Material 








(bulk tilt system) 


uj (cm ^) 


fc-poiiit 


distort ion- type 


LaAlOs 


46. 9i 


M 


a a c 




135i 


A 


a a c 




108i 


A 


r,- r,- 

CL a, c 


BiAlOs 


175i 


i 


^-polarization 




158i 


r 


a^y- polarization 




234i 


M 






211i 


It 


a, a c 




247i 


A 


a a c 




232i 






YAIO3 


96i 


r 


polarization 




248i 


M 


a'>a<'c+ 




229i 


R 


a+Q+c" 




291i 


A 


a°a°c- 




282i 


A 


a~ a~ 


LaGaOs 


mi 


M 


a'>a°c+ 


{a~ a~ c^) 


165i 


R 


a+a+c° 




215i 


A 


a°a°c~ 




209i 


A 


a~ a~ 


YGaOg 


180i 


r 


2-polarization 


(a-a-a+f] 


144i 


r 


a:y-polarization 




297i 


M 


a'>a°c+ 




286i 


R 


a+a+c° 




327i 


A 


a°a°c- 




322i 


A 


a~ a~ c'^ 



^ When confined to a perovskite manifold of structures. 



rion, indicating that a trilinear coupling is symmetry 
forbidden. 

(2) LaGaOa/YGaOa. — To mitigate this issue, we sub- 
stitute YAIO3 with LaGaOs to recover the trilinear 
coupling enabled by the chemical criterion. Since 
we find strong a'^a^c^ and a~a~c^ rotations present 
in both compounds (Table |l]) , we anticipate that 
the energetic criterion will be strongly satisfied — 
rotation-driven ferroelectricity should be active in 
the LaGaOa/YGaOa superlattice. 

(3) LaAlOa/YAlOa. — Starting again from (1), we will 
satisfy the chemical criterion by substitution of 
YGaOs with LaAlOa. In contrast to case (2), the 
energetic requirement in this all-aluminate super- 
lattice is weakly satisfied because LaAlOa displays 
the rhombohedral tilt pattern in its ground state 
rather than the targeted orthorhombic one; this 
is reflected in the substantially harder a°a°c+ and 
a~a~c'^ mode frequencies. 



(4) LaAlOa/BiAlOs. — Finally, rotation-induced ferro- 
electricity should be deactivated by substituting 
YAIO3 with BiAlOa because neither compound pos- 
sesses the orthorhombic Pnma tilt pattern — a direct 
violation of the energetic criterion despite satisfy- 
ing the chemical restriction. Proper ferroelectricity 
could emerge, however, from the strong polar insta- 
bility present in BiAlOa. 

Results 

We find that all superlattice compositions with A /A' 
cation ordering have polar crystal structures consist- 
ing of cation displacements and octahedral rotations 
(Fig. S3). The calculated ferroelectric polarizations 
shown in Table |ll] are comparable to conventional fer- 
roelectrics, e.g. the prototypical perovskite BaTiOa ex- 
hibits a ~33 /j,C cm~^ polarization in is ground state 
structure.!^ Consistent with our materials guidelines, the 
B/B' cation ordered superlattice is centrosymmetric. For 
the polar structures, however, the origin for ferroelectric- 
ity remains to be identified. The necessary condition to 
realize rotation-induced ferroelectricity is an interaction 
between different rotation distortions that produces a net 
energy gain to the polar ground state from their mutual 
coexistence. We determine this interaction by analyzing 
the energy surface around the paraelectric parent struc- 
ture in terms of the unstable lattice modes that produce 
the octahedral rotations present in the ground state struc- 
tures. We show from this analysis that complete control 
over hybrid improper ferroelectricity is gained through 
application of our simple design criteria. 

(1) YAlOa/YGaOa. — The centrosymmetric ground 
state structure exhibits anti-parallel cation displacements 
and the centric a~a~c~ tilt pattern. Here, adjacent oc- 
tahedra rotate out-of-phase in all directions and the xy- 
rotation angle magnitude is modulated from Ga-layer 
to Al-layer along the z-direction (Fig. S4). This octa- 
hedral motif results from a combination of the a°a°c~ 
and a~a~c^ unstable zone-boundary instabilities of the 
paraelectric superlattice. However, the combination of 
a'^a'^c^ (or the symmetry equivalent but energetically 
unique a^a'^c"'' rotation) with the a~a~c'^ mode in the 
B/B' superlattices produces only centrosymmetric struc- 
tures (Fig. S5) and therefore prohibits the octahedral 



TABLE IL Conditions for hybrid improper ferroelectricity 
(HIF) in 1/1 perovskite superlattices. For the octahedral 
rotations to induce a ferroelectric polarization in the cation 
ordered superlattices both the chemical {A/ A' ordering) and 
energetic (orthorhombic a~ a~ tilt pattern) criteria must be 
satisfied. 

A/A' polarization 
ordering? tilt pattern? (/iC cm"^) HIF? 

(1) YAlOs/YGaOa no no no 

(2) LaGaOs/YGaOs yes yes, strong 11.8 yes 

(3) LaAlOa/YAlOs yes yes, weak 9.94 yos 

(4) LaAlOa/BiAlOs yes no 9.62 no 
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a LaGaO:i/YGaO;, b LaAlOa/YAlO;, c LaAlOa/BiAlOa 
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FIG. 4. Energetic signatures for rotation-induced ferroelectricity. The calculated zero kelvin two-dimensional energy surface 
contours for each cation-ordered perovskite with respect to the primary Qi : and Q2 : 1 a c rotations centric octahedral 

rotation modes present in the polar ground state structures. Filled symbols indicate the positions of energy minima in the 
potential energy landscape. Equivalently colored symbols represent structures with identical ferroelectric polarization directions. 
Circles and squares differentiate anti-phase (domain) structures, which have identical polarization directions, but possess 
different non-polar structural distortion directions. In (a) and (b) the two and a a c rotations combine to produce the 

orthorhombic a~ a~ tilt pattern: (a) The strong rotational modes in the LaGaOs/YGaOs superlattice produce four polar 
crystal structures. Each Pmc2\ structure is symmetry related to the others as either ferroelectric twin structures, which differ in 
their polarization direction ([110]- or [110]-type), or as antiphase domains, which differ in the relative sign of the GaOs rotations 
while the polarization direction remains fixed, (b) Despite the rotation being substantially weaker than the a a c tilt 

in the LaAlOa/YAlOa structure, the combination of two non-polar octahedral rotations produce four symmetry related polar 
Pmc2\ structures. The energy barrier separating the ferroelectric twins, however, is substantially reduced, (c) In the cation 
ordered LaAlOs/BiAlOa superlattice this energy barrier collapses to zero: the aPaPc^ and a~ a~ cP strongly compete with each 
other to produce two symmetry equivalent non-polar (green triangles) Pmma structures, indicating the calculated ferroelectric 
polarization is not induced by the octahedral rotations. 



rotations from inducing ferroelectricity. Consistent with 
this result, we find that these rotations also produce an 
identical non-polar structure in the Bi substituted ana- 
logue, BiAlOa/BiGaOs (Fig. S6). Our first-principles 
predicted centrosymmetric P2i/c structures confirm the 
symmetry-derived chemical criterion: Compositions with 
B / B' cation ordering prohibit the centric octahedra from 
producing a ferroelectric polarization. 

(2) LaGaOa/YGaOa. — The ground state gallate super- 
lattice exhibits the targeted orthorhombic a~a~c^ octa- 
hedral tilt pattern and possess polar, rather than fully 
compensated antiparallel, cation displacements. Adja- 
cent GaOe octahedra rotate coherently in-phase (with 
the same magnitude) about the axis perpendicular to the 
La/Y ordering and out-of-phase in the two orthogonal 
directions (Fig. |3]) , which result from the combination of 
the highly unstable aPaPc^ and a~a~c^ rotations present 
in the paraelectric LaGaOa and YGaOa phases (Table |l]) 
that are each independently energy lowering. 

Rotation-induced ferroelectricity requires the two oc- 
tahedral rotation modes present in the polar structure 
to dominate the energy landscape. To explore this, we 
map out the two-dimensional energy surface contours for 
the LaGaOa /YGaOs superlattice in terms of these modes 
(FigureEa): We find four symmetry equivalent energy min- 



ima (denoted by filled symbols) with noncentrosymmetric 
Pmc2i crystal structures, indicating that the a'^a'^c^ and 
a~a~c'^ rotations combine to produce the ferroelectric 
structure in the absence of cation displacements. If the 
two modes were strongly in competition with each other, 
producing an energy penalty by their mutual coexistence 
in the system, then only two minima would result. The 
location of these minima deep inside the quadrants and 
the large energy barrier (~39 meV/f.u.) separating the 
ferroelectric twins (at Qi > and < 0) reflects the 
strong favorable coupling between the largely unstable ro- 
tational modes possessed by both bulk constituents. This 
barrier is further enhanced by nearly 3X by substitution 
of Y with Bi (Fig. S7). For these reasons, we anticipate 
rotation-induced ferroelectricity in this composition to be 
robust to thermal effects at room temperature. 

To verify that the a^a"c+ octahedral tilt pattern in- 
duces the electric polarization (and the subsequent cation 
displacements), we remove all polar distortions from the 
Pmc2i ground state structure and calculate the electronic- 
only contribution to the total polarization using the Berry 
phase method with only the octahedral rotations frozen 
in.^"^ Any electric polarization now must originate solely 
from the change in charge density induced by the coop- 
erative and non polar oxygen displacements (octahedral 
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rotations). Here, we find the rotations induce a sizable 
electronic-only polarization of 0.65 fiC cm~^. The po- 
larized charge density acts as an effective electric field, 
driving the ions to displace in a polar fashion — this is the 
essence of hybrid improper ferroelectricity. These induced 
cation displacements contribute to the total electric polar- 
ization in the fully relaxed Pmc2i structure. As a result, 
the polarization P is intimately linked to the a~a~c'^ tilt 
pattern in LaGaOa/YGaOa. 

(3) LaAlOa/YAlOa. — The ground state polar struc- 
ture of the all-aluminate composition is identical to the 
previous all-gallate superlattice: it exhibits the targeted 
a~a~c+ tilt pattern and polar cation displacements (space 
group Pmc2i). As before, we find four symmetry equiva- 
lent energy minima resulting from the combination and 
cooperation of the a'^a*'c+ and a~a~c^ rotations coop- 
erate with one another (Fig. |4j3). After full structural 
relaxation of this rotationally-only distorted structure 
(corresponding to the minima in Fig. [Sjo), we find po- 
lar cation displacements occurring in response to the 
a~a~c+ tilt pattern through the symmetry allowed tri- 
linear Pz{QiQ2) coupling term absent in single phase 
non-polar perovskites. 

Although the net interaction between the two rotation 
modes lower the energy of the LaAlOs/YAlOa system, Q2 
is stronger than Qi. This is evident from inspection of the 
unstable phonon modes for the bulk constituents (Table |l| 
and is also discernible in Fig.[4]D, where the energy minima 
are located close to the Qi = boundary. This weak 
susceptibility to the aPaPc^ tilt system in paraelectric 
LaAlOs, leads to a small (~1 meV/f.u.) zero-kelvin energy 
barrier separating the two ferreoelectric twins. Thus 
the energetic balance between the two rotation modes 
and degree to which the energetic criterion is satisfied 
is equally as important as maintaining an A/ A' cation 
ordered superlattice (chemical criterion). It determines 
the stability of octahedral rotation-induced HIF: Because 
the aPaPc^ rotation instability in LaAlOs is weak, its 
presence in the ground state crystal structure is highly 
susceptible to finite temperature effects. If the four polar 
minima are thermally suppressed then two non-polar 
Pmma crystal structures would result, subverting the 
energetic criterion and nullifying the rotation-induced 
ferroelectricity mechanism. 

(4) LaAlOa/BiAlOa. — The ground state polar struc- 
ture exhibits the targeted a~a~c^ tilt pattern; however. 
Figure [4j; illustrates that in the absence of polar ionic 
displacements, the energy interaction between the two 
aPa^c^ and a~a~c^ rotations disfavors their mutual 
coexistence. Only two symmetry equivalent energy 
minima with centrosymmetric Pmma structures occur. 
In this extreme case, the energy cost of the coupling 
between rotations exceeds the gain from the individual 
modes. In other words, while the rotations Qi Q2 are 
of suitable symmetry to induce an electric polarization 
and are individually unstable (Table S4), our density 
functional calculations reveal that one mode dominates 
and suppresses the other. This happens in part due to the 



weak Qi instability in LaAlOa and strong anharmonic 
coupling between Qi and Q2 in bulk BiAlOa. Both 
of these features are reflected in the fact that neither 
lanthanum aluminate nor bismuth aluminate display the 
desired a~a~c+ tilt system in their ground states. The 
ferroelectric polarization that results in this system is 
due to a conventional soft polar mode; rotation-induced 
ferroelctricity is absent in the LaAlOa/BiAlOs system. 

Outlook and summary 

Finally, we describe important applications exploiting the 
trilinear coupling that enables centric octahedral rotation- 
induced ferroelectricity. This mechanism supports electric 
polarizations in more diverse chemistries. Unlike the con- 
ventional SOJT-mechanism that lifts inversion symmetry 
and produces electric polarizations through cooperative 
displacements of cations with dP electronic configurations 
(group 4, 5 and 6 transition metal ions) or cations with 
stereochemical active ns^ lone pair electrons, no such re- 
striction is imposed on the cations' valence in this form of 
hybrid-improper ferroelectricity. Cations with strong mag- 
netic interactions, open d- or /-shell configurations, which 
are incompatible with conventional ferroelectricity,'^ are 
able to both fully coexist and couple to sizable electric 
polarizations: Centric octahedral rotation-induced ferro- 
electricity provides a plausible route to achieving robust 
magnetoelectric multiferroics. 

It is interesting to conjecture that there are viable 
chemistries providing rotation-induced ferroelectrics with 
anomalously large electronic- only contributions to the 
total polarizations. In cases where the polarization con- 
tribution (and distortions) due to ionic displacements is 
also small, reversal of the electric polarization would re- 
quire minimal ionic motion. These materials would satisfy 
many of the robust cycling and ultra fast (femtosecond 
timescale) dynamical switching requirements needed for 
next-generation electronics. 

We have used first-principles density functional 
calculations combined with group theoretical studies 
to enumerate the crystal-chemistry design criteria 
required for the centric octahedral rotations pervasive 
in perovskite oxides to induce ferroelectric polarizations. 
We showed that interleaving two bulk perovskites to 
form an ordered and layered arrangement of A /A' 
cations [chemical composition {A, A')B20q\ produces 
a new trilinear free energy term coupling three lattice 
modes: two octahedral rotations that dominate the 
energy landscape, forming the orthorhombic a~ a~ c^ 
perovskite tilt system, and an electric polarization. For 
ordered perovskites where these conditions are satisfied, 
the trilinear term induces an electric polarization and 
hybrid improper ferroelectricity results. Because most 
octahedral rotations in perovskites freeze-in at elevated 
temperatures (>300 K), we argue that the trilinear 
coupling provides a robust route to realize synthetic ferro- 
electrics from constituents which are not ferroelectric in 
the bulk at room temperature. By leveraging systematic 
symmetry analysis with first-principles density functional 
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calculations, we illustrated an a priori materials selection 
strategy for designing synthetic ferroelectric crystals only 
from knowledge of the lattice dynamics of the constituent 
materials. 

Methods 

Our zero kelvin density functional calculations are per- 
formed within the local density approximation (LDA) 
using the Vienna ab initio Simulation Package (vASpj^^^ 
and the projector augmented wave (PAW) methocPl 
with the following valence electron configurations: 
5s^5p^5d^6s^ (La), 4524/5524^1 (Y), 3d^°4:S^4p^ (Ga), 
Ss^Spi (Al) and 2s^2p'^ (O). For the reference Pi/ mmm 
bulk perovskite structures we use a 5 x 5 x 5 Monkhorst- 
Pack fc-point meslpS' and a 500 eV plane wave cutoff. We 
then volume-optimize each structure within PA/ mmm 
symmetry. For our superlattice calculations we use a 
5x5x3 Monkhorst-Pack fc-point mesh and a 550 eV 
plane wave cutoff. We relax the ionic coordinates until 
the Hellmann-Feynman forces on the atoms are less than 
0.1 meV A-i. 

We obtain the phonon frequencies at high symmetry 
Brillouin zone points (Table [l| by calculating total en- 
ergies with respect to atomic displacements about the 
reference PA/ mmm structure in a 2x2x2 supercell. For 
these calculations we use a larger, 650 eV, planewave 
cutoff. In this frozen-phonon method, a series of small 
(symmetry inequivalent) atomic displacements are im- 



posed along different Cartesian directions. We calculate 
the dynamical matrix from the Hellman-Feynman forces 
and diagonalize the dynamical matrix to obtain the atomic 
displacement patterns (eigenvectors) and phonon mode 
frequencies (eigenvalues) . 

We obtain the ground state atomic structures for all 
cation ordered perovskite compositions by systematically 
"freezing-in" linear combinations of the unstable phonon 
modes (Table |l]) into the paraelectric cation-ordered 
perovskite reference structures (Sec. SII.B) and then 
performing full structural relaxations from first-principles. 
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